CanopyCover Function

public function CanopyCover(dc, den) result(cc)

compute canopy cover (0-1)

References:

Collalti, Alessio & Perugini, Lucia & Santini, Monia & Chiti, Tommaso & Nolè, Angelo & Matteucci, Giorgio & Valentini, Riccardo, 2014. A process-based model to simulate growth in forests with complex structure: Evaluation and use of 3D-CMCC Forest Ecosystem Model in a deciduous forest in Central Italy, Ecological Modelling, 272(C), 362-378.

Arguments

Type IntentOptional Attributes Name
real(kind=float), intent(in) :: dc

crown diameter (m)
real(kind=float), intent(in) :: den

tree density (trees/ha)

Return Value real(kind=float)

returned result


Variables

Type Visibility Attributes Name Initial
real(kind=float), public :: dbhdc

actual ratio between dbh and crown diameters (m/cm)

Source Code

FUNCTION  CanopyCover &
!
(dc, den) &
!
RESULT (cc)

IMPLICIT NONE

!arguments with intent(in):
REAL (KIND = float), INTENT(IN) :: dc !! crown diameter (m)
REAL (KIND = float), INTENT(IN) :: den !! tree density (trees/ha)

!local declarations:
REAL (KIND = float) :: cc !! returned result
REAL (KIND = float) :: dbhdc !! actual ratio between dbh and crown diameters (m/cm)

!-----------------------------------end of declarations------------------------


cc = Pi * dc ** 2. / 4. * den / hectare

!check boundary
IF ( cc > 1. ) THEN
    cc = 1.
END IF

IF ( cc < 0.) THEN
    cc = 0.
END IF


RETURN
END FUNCTION CanopyCover